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ABSTRACT 



The developing flow in the hydrodynamic entrance region of an equi- 
lateral triangular channel was investigated 0 Using a fully developed 
velocity profilfi from Knudsen and Katz (3) as a boundary condition s the 
equations were solved numerically employing the method of Chorin (1) * 

The resulting velocity profiles and pressure drop were employed 
in calculating the local friction factor in the hydrodynamic entrance 
region* 
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I® INTRODUCTION 



The flow characteristics of non-circular channels have become 
of increasing importance in the analysis of compact heat exchangers® 
Nevertheless, heat transfer analyses of the triangular entrance region 
have been limited to approximations which postulate fully established 
velocity profiles (Graetz Approximation)® Such solutions have been 
obtained by Lu and Miller (4) and Tao (8) for equilateral triangles® 
Using Sparrow’s (7) results for the fully established velocity 
profiles in an isosceles triangular duct, McComas (5) has estimated 
values of entrance lengths and friction factors for such flows® 

In order to complete the analysis of heat transfer in the entrance 
region of triangular ducts, it is necessary to obtain a hydrodynamic 
solution for the developing velocity profiles in that region® It is 
therefore the purpose of the present work to provide a solution of 
the equations of motion in the hydrodynamic entrance region of an 
equilateral triangular ducto 
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lie ANALYSIS 



In an attempt to find a convergent solution to the equations of 
motion* several explicit methods were tried without success » Subse- 
quently a relaxation method employing an artificial density change* 
suggested by Chorin (1)* was adopted© This method essentially alters 
the solution of the equations of motion from an initial value problem 
to a boundary value problem* necessitating the use of entrance length 
and the fully developed velocity profile as boundary conditions© 

Ac Governing Equations 

The governing equations for the flow of an incompressible viscous 
fluid (Navier Stokes Equations) may be written: 
x - Momentum: 



^ + U ^ * V ^ + ur ^ - - i *=P U 

^ cJx ^ ^ Ux 



( 1 ) 



y - Momentum: 

u + v ^ ur = _ j + -o u- 
^ f ^ 



z - Momentum: 

4- U ^ 4 \r ^ -v ur ^ - - 1 V z uT 

•it ^ ^ 



Continuity l 



4 * + 





( 2 ) 



( 3 ) 



( 4 ) 



12 



Introducing the dimensionless variables. 



u- 


u. 

U a 


V- 


v~ 

Uo 


w* 


X s 


X 

d 


Y a 


a 
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7 = 



P = -p 



?l)Uo 



ur 

Uo 



t =t 



T> 

d* 



equations 1-4 become: 






it 



, + R 



iV 



it' 



: + r 









+ R fu ^ + v +W 






j- 

]- 



it? 

3Y 



v v*V 



^ p + V 2 w 

52 



*V iW _ n 

JX + 5y * 57 "° 



( 5 ) 



( 6 ) 

(7) 

(8) 

(9) 



B» Boundary Conditions 

A condition of uniform velocity is assumed at the entrance and 
no-slip conditions are imposed at the walls; thus: 

1. at the entrance:* 

tr(o,Y,?,t) » 1 
V (0,Y,7,t) - 0 
W(0,Y,?.t) • 0 

2* at the walls i 

U(t)*V(i)'W <L\°0 



*The prime has been dropped and dimensionless variables are 
understood^ 
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Co Artificial Density 



Because the problem is one of incompressible flow, no density 
change can take place 0 One can introduce an auxiliary continuity 

equation! ^ _ 

^ * 



(9a) 



in which changes of the pseudo density given by 



P= F/S (10) 

can be used to test for convergence at each iteration process » 

Thus, we have introduced a set of auxiliary parameters, ^ and £ , 

which are analogous to an artificial density and an artificial com- 
pressibility 0 

By introducing the artificial equation of state (equation 10) 
and using t* as an auxiliary parameter analogous to time in compressible 
flow* the initial value problem can be altered to a boundary value 
problem© Such a change, however, requires the use of a fully estab- 
lished solution as an additional boundary condition An exact solution 
for the steady, fully established profile has been obtained by Knudsen 
and Katz (3) : 

IUY.2) - -2fR [$(?>♦ 3?T*)-4] (u) 

where? 7-Y~jy^and2 ? ~ 3d ° The friction factor f for fully 

developed flow is given by Kays (2) as 13©33 0 This number serves to 
determine the value of the fully established pressure gradient© 

Do Finite Difference Equations 

Equations 6-9 were put into finite difference form using central 
differences, except for the second derivatives which were differenced 
using the Dufort-Frankel pattern 0 This latter expression is of the forma 
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( 12 ) 



tl 

ST 






xg ^ u cx-o » if (x-) - u -u 

rt aX z 
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where U N (X+) denotes U(N,X+AX,Y,Z), uN+1 is U(N+1,X,Y,Z) and N 
denotes the Nth relaxation e The resulting equations were solved for 
the N+l terms and the resulting computing equations are of the form: 

u~' - - $ )] « {-& i [*“«-> - ^-0 

- Rfj [uS*,U «.| - - R$ [uS-)Vw*i - Uw-)V«-)] 

“ R jue*)W&+> - UJa-iWF* -->J -*• [u'w* ’ " kj J 



v“*'= W(**5Jk&*&)} x 

" [v wb>)-Yix-)US;o] - Rfy [V^ +) VCf^“ 

- Ra fv&oW(f^-v&oWi-)J+ a 4 * [v<^>+vS-,- v“" J 

4- [v^v&o-v*-’]} 



(13) 



(14) 



w* 1 ’ [Vh + iT*^"s‘)) x t'&Hf®**- 

-R2ljy«*>Uw>-WS-)l)«-0 - Riv(W«*)Viw-'Wo'->V«-)] 
~ R ^ jWa*)Wte+) -'VJte-)Wci-'] + fW (*■*> 4 'W (x- ) - J 

+ |y /r ^ + > ^'W"<Y-> -~W' J j 4- 2^ jW(2^)4--W&-,- W >J 'jl 



^ 



(J-i 



ju5*) - Ik*-) ] 



— ^ |Y(x+) -V(v->j - -"Wt? >J 



(16) 
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E 0 Grid Size and Stability 



The geometry of an equilateral triangle presents a problem in 
selecting appropriate values of AY and AZ* The requirement that the 
network intersections fall on a boundary was considered essential for 
numerical calculation* dictating the relationship between aY and AZ 
such that aY - V3 * A? /Z. . 

The value of AX, while arbitrary, is constrained by the fact 
that the number of X-wise stations must be chosen such that the number 
of stations times the value of AX be equal to the predicted entrance 
length, or AX x I=X e , where I is some integer© Finally, because of 
the limited storage available in a computer, the number of grid points 
is limited© For the purposes of this investigation, the values adopted 
were aX=©C0159, aY= s 0©17321, and aZ= 0©3© This grid resulted in 

25 stations in the X-direction, 11 in the Y and 6 in the Z f as shown 
in Figure 1 0 

Having selected the grid size, the stability of equations 13-16 
may be determined using the results of Chorin (1), who points out that 
the introduction of an artificial density results in an artificial 
speed of sound of the forms 

C2 " /§ (17, 

A Mach number relative to the speed of sound is given by, 



For convergence it is necessary that the Mach number (M) be less than 
oneo With this condition satisfied, the system of equations is stable 
when! 




(18) 




( 19 ) 
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In the present work these values were taken to be >1=00975, §=9*38025, 

At^O* 00035, and R=l*0 o 
F* Friction Factor 



The wall shear stress may be related to the local pressure gradient, 
which is of importance to engineers in predicting pressure drop* The 
average stress at any section is given by the Newtonian shearing law, 



Using equation 20, a dimensionless friction factor can be defined, 



in which f is the usual flow friction factor* 

G* Numerical Solution 

In the present analysis, using the values of AX, aY, aZ, At 
and § given above, equations 13-16 were found to converge to a steady 
solution after 2,000 iterations* While Chorin points out that initial 
values of variables may be set at zero throughout, the convergence can 
be hastened by using a good "guess" at the values to be calculated* 

In the present case the fully established velocity values were used 
throughout as an initial assumption* 

A copy of the program written in Fortran IV language is given in 
the Appendix* 




( 20 ) 



where 




represents the average velocity gradient evaluated at 



the wall* 




( 21 ) 
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IIIo RESULTS 



The dimensionless centerplane velocity profiles are given in 
Figure 2 for five values of X/X e (0*2, 0©4 f 0 o 6 c 0©8 f and 1©0) 0 Figures 
3 through 7 show the centerplane velocity profiles and velocity contours 
for X/X e =0©2 9 0©4 t 0©8 P and l o 0 respectively © The secondary ve- 

locities 9 V and W f are shown in Figure 8 at a station about midway down 
the entrance region at X/X e »0©6 6 The value of friction factor at each 
x-wise station is plotted in Figure 9 in the form of f/f e versus X/X e> 
where the e subscript denotes the fully established values© 

The dashed line on the friction factor curve was computed from 
the local skin friction predicted from a boundary layer analysis© The 
value of skin friction must of course be unbounded at X=0 because of 
the discontinuity in velocity at the lip of the duct© For large values 
of X/X €e f/f e must approach unity asymptotically as the results indi- 
cate© 



f 
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IV. CONCLUSIONS 



From the results of this study * it may be concluded that: 

1. The relaxation procedure used in this analysis was found to 
be convergent subject to certain stability parameters as Chorin (1) 
points out tt The results* however* remain unchecked experimentally. 

2. The use of the fully established values as exit boundary 
conditions forces the relaxation values to converge on these values 
in the limit as the number of iterations increase. Therefore* the 
values of the computed velocity profiles can be only as accurate as 
the exit boundary conditions. Due to the lack of analytic and experi- 
mental results in the area* there was little choice in using the 
results of Knudsen and Katz (3) for the fully established velocity 
profile and Kays (5) for the fully established pressure gradient. 

3. In general* the grid sizes used should be reduced until the 
number of points reaches the limit of computer storage. This serves 
to maximize the accuracy of the results consistent with the accuracy 
of the fully established profiles used 0 

4. It seemed intuitively obvious when first using Chorin’ s method 
that because it represents a convergent numerical solution to the 
finite difference form of the Navier Stokes and continuity equations* 
using an entrance length greater than a predicted value would lead to 
the fully established profile upstream of the boundary condition. 

This* however* was not the case 0 Thus* the necessity of using the 
predicted value of X e reduces the generality of this method in the 
solution of physical problems 0 The use of a Reynolds number of unity 
may have been the cause of this* however it would seem evident that a 
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complete analysis of the equations must be made numerically to find 
the reason for this limitation 

5 0 The Reynolds number dependency of the governing equations can 
be extracted by using non-dimensional space parameters of the form 
X=xR/d 0 An analysis of the resulting equations would seem beneficial 
as the results would be completely general and would represent "uni- 
versal" profilesc 
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FIGURE 2. 

DIMENSIONLESS VELOCITY PROFILES 
FOR VALUES OF X/X # . 
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FIGURE 3. 

VELOCITY PROFILE AND CONTOURS 
FOR X/X # =0.2 

24 





FIGURE 4. 

VELOCITY PROFILE AND CONTOURS 
FOR X/X # *0.4 



25 





VELOCITY PROFILE AND CONTOURS 
FOR X/X # s 0.5 
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FIGURE 6. 

VELOCITY PROFILE AND CONTOURS 
FOR X/X =0.8 
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FIGURE 7. 

VELOCITY PROFILE AND CONTOURS 
FOR X/X # =I.O 
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FIGURE 8. 

SECOND VELOCITY PROFILES 
FOR X/X e = 0.6 
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FIGURE 9. 

FRICTION FACTOR RATIO vs. X/X e 
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